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Abstract 

We investigate generalizations of the Charlier and the Meixner polynomials on 
the lattice N and on the shifted lattice N + 1 — /3. We combine both lattices to 
obtain the bi-lattice N U (N + 1 — /3) and show that the orthogonal polynomials 
on this bi-lattice have recurrence coefficients which satisfy a non-linear system of 
recurrence equations, which we can identify as a limiting case of an (asymmetric) 
discrete Painleve equation. 
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1 Introduction 

The classical orthogonal polynomials (Jacobi, Laguerre and Hermite polynomials) are 
orthogonal polynomials on (an interval) of the real line with a weight function that satisfies 
a first order differential equation (the so-called Pearson equation) of the form 

[a{x)w{x)]' = r{x)w{x), 

where a is a polynomial of degree at most 2 and r is a polynomial of degree 1. This 
Pearson equation allows to find many useful properties of these polynomials. There are 
other families of orthogonal polynomials in the Askey table [5J which live on a linear lattice 
rather than on an interval. These are the Charlier polynomials (O §9.14 on p. 247], [3l 
Chapter VI, §1]) 

OO h. 

fc=0 
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which are orthogonal on the lattice N with respect to the Poisson distribution, and the 
Meixner polynomials ([S', §9.10 on p. 234], [3, Chapter VI, §3] who calls them Meixner 
polynomials of the first kind) 

J2 Mn{k; /3, c)Mmik; /3, c)^^ = (^y^^Y^^*^-,- /3 > 0, < c < 1, 



k=0 



which are orthogonal on N with respect to the negative binomial (or Pascal) distribution. 
The special case /3 = — and c = p/{p — 1), with N & N and < p < 1 are the 
Krawtchouk polynomials which are orthogonal on the integers {0, 1, 2, ... , A^} with respect 
to the binomial distribution. The Hahn polynomials are another family of polynomials 
which are orthogonal on the integers {0, 1, 2, ... , A^}. Instead of the differential operator 
it is better for these lattice polynomials to work with the difference operators 

V/(x) = /(x) — /(x — 1), (backward difference) 

and 

A/(x) = f{x + 1) — /(x), (forward difference). 
The weight = /k\ for Charlier polynomials now satisfies the Pearson equation 



Vw(x) = ^1 j w{x) 



if we define the weight function w{x) = a^'/r(x + 1), so that Wk = w{k). This is of the 
form V[a{x)w{x)] = t{x)w{x) with a = 1 and r a polynomial of degree 1. The weight 
for Meixner polynomials Wk = {(3)kC^/k\ can be written as Wk = w{k) using the weight 
function w{x) = T{f3 + x)c^/(r(/3)r(x + 1)) and this satisfies the Pearson equation 



+ x)w{x)] = (^(3 + X j w(x). 



so that cr(x) = (3 + X and r is a polynomial of degree 1. This Pearson equation (with 
a difference operator) allows to find many properties of these orthogonal polynomials on 
the lattice M. In particular one can find the three term recurrence relation for these 
polynomials explicitly: 

— xC„(x; a) = aC„+i(x; a) — {n + a)C„(x; a) + raC„_i(x; a) 

and 

(c - l)xM„(x; /3, c) = c{n + /3)M„+i(x; /3, c) - [n + {n + /3)c]M„(x; /3, c) + nM„_i(x; /3, c). 

Often it is more convenient to work with the monic polynomials P„(x) = x" + ■ ■ ■ , for 
which the recurrence relation is 

XP„(X) = Pn+l{x) + hnPn{x) + a^P„_l(x), 

or with the orthonormal polynomials Pn{x) = In^^ + ■ ■ ■ ? with 



oo 

In 



k=0 
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for which the recurrence relation is 



xpn{x) = a„+ip„+i(x) + bnPn{x) + a„jo„_i(x). 



The recurrence coefficients for Charher polynomials then can be seen to be = na and 
bn = n + a for n G N, and for Meixner polynomials one has 



In this paper we will investigate generalizations of the Charlier and the Meixner poly- 
nomials. In Section 2.1 we will use the weight function w{x) = r(^)a^/(r(/3 + x)r(x + l)) 
which gives the weights Wk = w{k) = / on the lattice N. In Theorem 12. II we will 
find a system of non-linear recurrence equations for the recurrence coefficients a\ and 6„ 
and we can identify this system as a limiting case of a discrete Painleve IV equation. In 
Section 2.2 we will use this weight on the shifted lattice N + 1 — /3 and in Theorem 12.21 we 
will show that the recurrence coefficients satisfy the same non-linear system of recurrence 
equations, but with a different initial condition for b^. In Section 2.3 we combine both 
lattices to obtain the bi-lattice N U (N -|- 1 — The orthogonality measure is now a 
linear combination of the measures on N and N -|- 1 — /3 and surprisingly, Theorem 12.31 
shows that the recurrence coefficients again satisfy the same non-linear system of recur- 
rence equations, but with a different initial condition for b^. In all these cases, the initial 
conditions for bo are given by ratios of modified Bessel functions. In Section 3 we deal 
with a similar generalization of the Meixner weight and we investigate the weight function 
w{x) = r(/3)r(7 -|- x)a^/(r(7)r(/3 + x)T{x + 1)). In Section 3.1 we use this weight on 
the lattice N where it gives the weights Wk = w{k) = {pf)ka^ / {{[3)kk\). In Theorem 13.11 
we will show that the recurrence coefficients of the corresponding orthogonal polynomials 
satisfy another non-linear system of recurrence equations, which we can identify as a lim- 
iting case of an asymmetric discrete Painleve IV equation. In Section 3.2 we use the same 
weight function but now on the shifted lattice N+ 1 — /3. The recurrence coefficients of the 
corresponding orthogonal polynomials satisfy the same non-linear system of recurrence 
equations, but with a different initial value for b^. Finally, in Section 3.3 we study the 
orthogonal polynomials on the bi-lattice N U (N + 1 — /3) with respect to a linear combi- 
nation of the measures on N and N -|- 1 — /3, and again the recurrence coefficients satisfy 
the same system of non-linear recurrence equations, but with a different initial value for 
6o- The initial conditions in Section 3 are all given by ratios of confluent hypergeometric 
functions. 

We hope this paper is of interest to people studying orthogonal polynomials because 
we introduce some new families of orthogonal polynomials which seem to have nice prop- 
erties. The asymptotic behavior of the recurrence coefficients may be of interest since 
in many cases one can observe an oscillating behavior, reminiscent of orthogonal polyno- 
mials on two intervals. The paper should also be of interest to people studying discrete 
Painleve equations since we are giving two systems of discrete Painleve equations and a 
one-parameter family of initial conditions for which the solution turns out to be in terms 
of recurrence coefficients of orthogonal polynomials on the real line. This means that 
these particular initial values do not lead to singularities and in fact have nice positivity 
properties, such as a„ > 0. The fact that the polynomials are orthogonal on a (bi-)lattice 
also gives some properties for partial sums of the bn (sums of the zeros). The cases where 




'n 




1 - 



C 
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the bi-latticc reduces to a single lattice may well be corresponding to unique solutions 
of the discrete Painleve equations with a prescribed asymptotic behavior for the 6„ as 
n — >■ oo. 



2 A generalized Charlier weight 

2.1 The lattice N 

We consider the discrete weights 



*^N={0,1,2.3....} 
with a > 0. We can write these in terms of the weight function 

which is now a function of a; G C vanishing whenever a; is a pole of r(x + 1) or a pole 
of r(/3 + x), i.e., for x = —1, —2, —3, . . . and x = — /3, — /3 — 1, — /3 — 2, . . .. This weight 
satisfies the Pearson equation 

Vw(x) w(x) - w(x (2.2) 

a 

With this weight we can introduce the inner product 

oo 

which has positive weights for every /3 > 0. We denote the corresponding monic orthogonal 
polynomials by Pn{x; a, (3): 

oo ^ 

Y^P^{k;a,P)PUk;a,/3)-^ = 0, n ^ m. (2.3) 

fc=0 

Our main interest is to find the recurrence coefficients in the three term recurrence relation 

xPn{x;a,/3) = Pn+i{x;a, /3) + bnPn{x;a, /3) + alPn-i{x; a, /3), (2.4) 

with initial conditions P_i = and Pq = 1. We will obtain recursive relations for these 
recurrence coefficients. We sometimes just write Pn{x) for the polynomials Pn{x; a, /3) to 
simplify the notation. Observe that for a — f3c and ^ — > oo we get the weights Wk — /k\ 
so that 

lim Pn{x;/3c,/3) = C'„(x;c) 

are the monic Charlier polynomials. 

The Pearson equation and summation by parts gives the following structure relation 
for the orthogonal polynomials 
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Lemma 1. The monic orthogonal polynomials given by satisfy the relation 

APn{x) := Pn{x + 1) - Pn{x) = nPn-i{x) + 5„P„_2(x), (2.5) 
for some sequence Bn of real numbers. 

Proof. We can always write the polynomial Pn{x + 1) — Pn{x) (which is of degree n — 1) 
in a Fourier series using the orthogonal polynomials and hence 



n-l 
k=0 



Pn{x+l)-Pn{x) = Y,Ak,nPk{ 

where the Fourier coefficients are given by 

{Pk,Pk) ■ 

Recall summation by parts: if 6_i = then 

CO oo 
A:=0 k=0 

whenever a and b are in £2- If we apply this, then 



{Pk,Pk)Ak,n = J2^AP-{j))PkU) 



Wj 



j=Q 



j=0 
00 



-J2Pnij)WjVPk{j) -J^PnimU - 1)V 

j=0 j=0 



Wj. 



The first sum on the right is zero by orthogonality since VPfc is a polynomial of degree 
k — 1 < n. For the second sum we use the Pearson equation (12.2p to find 



(Pfc, Pk)Ak,n = —y2 Pn{3)Pk{3 -l){a- 3(13 - 1) 



j=0 



This sum is zero by orthogonality whenever k + 2 < n. Hence only An-i^n and An-2,n can 
be non-zero. By comparing the leading coefficients, we see that An-i,n = n. If we call 
An-2,n = Bn, then the required formula follows. □ 

Some of these coefficients are useful in later computations. If we define 

9 {P'n, Pn) y 
7n 

then it is not difficult to find from the recurrence relation (12. 4 p (taking the inner product 
with Pn-i) the well-known relation 

2 In-l 



7n 



Furthermore, with the notation in the proof of the lemma 

1 = ^,1 = E - - 1) - fH. 

(Jj 

j=0 

By orthogonahty we have 

oo 
3=0 

If we use the recurrence relation then 

oo oo 

E jPi(j>, = J2[P2{j) + &iPi(j) + alPo{j)H = «?/7o- 

j=0 j=0 

If we use the recurrence relation twice, then 

oo oo 

Y^fPiUH = J2^h{P,{j) + hP^{j) + alPo{j)) + al{P,{j) + boPom 

j=0 j=Q 

= al{bi + bo)/jl 
Hence, combining these results, we have 

l = -(&i + &o + /3-l). (2.6) 

a 

In a similar way we can compute B2: 

2 

B2 = Ao,2 = 5^P2(j)(« - - 1) -f)wr 

By orthogonality 



i=o j=o 
and using the recurrence relation we get 



J2fP2{j)w, = Y,j{P,{j) + h2P2{3) + alP,{j)) 

2 2/2 

a2«i/7o 



3=0 3=0 

,2„2 /.,2 



SO that 

P2 = (2.7) 
a 

We now are able to derive the recursive relations for the recurrence coefficients. 

Theorem 2.1. The recurrence coefficients for the orthogonal polynomials defined by ^2.3\) 
satisfy 

an 

bn + bn^i-n + (3 = — , (2.8) 
(a^+i - a)(a^ - a) = a(&„ - n)(6„ - + /3 - 1), (2.9) 
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with initial conditions 

= 0, 60 



where 1^ is the modified Bessel function 



k=0 

Proof. On one hand we have the three-term recurrence relation (12.4p 

kPM = P^+,{k) + hMk) + alP„,^^{k), 

and on the other hand we have the structure relation (12. 5p . The compatibility between 
these two equations will result in the desired equations. Take the forward difference of 
the three-term recurrence relation to find 

(fc + 1)A(P„(A;)) + Pn{k) = APn+,{k) + bnAPnik) + alAPn-i{k). 

Then use the structure relation (12. 5p to find 

{k + 1) {nPn-l{k) + BnPn-2{k)) + P„,(fc) 

= {n + l)P^{k) + 5„+iP„_i(A;) + hn {nPn-i{k) + 5„P„-2(A;)) 

+ al {{n - l)Pn-2{k) + P„-iP„-3(A;)) • 

Finally, use the three-term recurrence relation to work out kPn-i{k) and kPn-2{k) on the 
left hand side. This gives a linear identity involving Pn{k), P„_i(/c), P„_2(fc), P„_3(/c), and 
since these four polynomials (in the variable k) are linearly independent, the coefficients 
should all be zero. This gives the four equations 

n + 1 = n + 1, (2.10) 

n + nbn-i + Bn = P„+i + (2.11) 

Bn + nal_^ + = Br^bn + {n - l)al, (2.12) 

Bna^_2 = a^Bn-i- (2.13) 

Clearly (12.101) is always satisfied and (I2.13P readily gives 

Bn _ Bn-l 

SO that Bn = a^a^_iB2/ (alal), and hence (12. 7p gives 

Bn = ^^I^, n>2. (2.14) 
a 

Use this in (1211]) to find 

na {bn - bn-i - 1) = al{al_^ - a^+J, (2.15) 

and (I2.12P becomes 

'-(b. - b._. - n = - 



Summing the latter starting from n = 2 gives 



1 N 1 /, 7 \ ''^ 1 

-[On + bn-i -n + 1) (6i + bo) = ^ 2- 

a a af 

Now use (EJ]) to find ( M . 

If we use (12 .Sp in fl2.15p and put 6„ = n + then we find 

(4 + 4-1 + A; + /3 - l)(4-i - 4) = flfc+i - al_^. 
Summing from k = 1 to n gives (use the telescoping property and summation by parts) 

n-l 

- dl + dl + ^dk - ndn - {(3 - 1)(4 - 4) = + al - aj, (2.16) 

fc=0 

where we used the initial condition = 0. On the other hand, (I2.15P is equivalent with 

a/c(4-i — 4) = o-fcO-fc+i — ctfc-iOfc- 
Summing from A; = 1 to n now gives 

n-l 

a 4 — a^4 = «n«n+l- (2-17) 

fc=0 

If we use fl^rrp in fl^T^ then we find 

-dl + dl + - (/3 - 1)(4 - 4) = a^+i + - a?. (2.18) 

The initial values c/q = feo and are given by 



mi 2 ^2 I TTT'l 



bo = — , = — 

mo mo \mo 



where mj are the moments 



oo 

m^ 

k=0 

A simple calculation gives 



where and Ip-i are modified Bessel functions. This gives 



_ /- //3(2v/^) 



The Pearson equation (12. 2p gives 



A:=0 fc=0 fc=0 
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which readily gives m2 = arriQ — {13 — l)mi. Observe that this gives 



c/q + (/3 — l)do + al = a 
which simphfies the recurrence relation f l2.18p to 

2 2 

2 I 2 I t2 .in i \ r n 

a„+i + a„ + ct„ V \P - l)"n - a = 0, 

which is equivalent with (12. 9p . □ 

The case /3 = 1 was already considered in [9] (see also [H §4.2]). Equation (12. 9 p then 
simplifies to 

{al^^-a){al-a) = {bn-n)^ (2.20) 

so that — a and a^+i — a have the same sign. The boundary condition Og = thus 
implies that — a < for all n, and we can write — a = — ac^, for some new sequence 
(cn)nGN- Then = a(l — c^) so that < 1 for n > 1. This still leaves two choices for the 
sign of Cn- Taking square roots in (I2.20p gives bn = n + aCnCn+i, where we choose Cq = 1 
and we recursively take the sign of c„+i equal to the sign of (6„ — n)/cn- Inserting these 
formulas for and hn into (12. 8p with (3 = 1 gives a(c„+i + Cn-i) = n/{l — c^), which is 
the discrete Painleve II equation ([1], O Appendix A.l]). 

The situation is different for /9 7^ 1, since (12. 9 p now is a quadratic equation in 6„ and 
it is not a priori clear which of the two roots one should choose. Equations (I2.8p - (l2.9p 
are a limiting case of a discrete Painleve IV equation: take the first dPiv in [8, p. 723], 
or equivalently the second discrete Painleve equation for in p. 297] 



Vn + Vn-l 



yl-B 

C , Cn + C 



1 + DXn 1 + Xn/D 

where = + n5 and Cn = Zn — 8/2. If we put x„ = iXn/ \faB and let 5 —> 00 then 
this gives for the first equation 

Xn+lXn = a ((?/„ - Znf - A) , 



and if we take iD = y^B/a and let 5 —i- 00 then we also get 

Vn + Vn-l = , +Cn + C. 

1 + X„/a 

The parameters A = {(3 — 1)^/4, C = —f3/2, Zn = n — {(3 — l)/2 then give the discrete 
equations (I2.8p -( l2l9l) for X„ = — a and ?/„ = 6„. 

2.2 The shifted lattice N + 1 - /3 

We can consider the weight w in (12. ip also on the shifted lattice N + 1 — /3 = {1 ~ (3,2 — 
/3, 3 — /3, . . .}, where 

9 



The weights {vk)ke'N are therefore of the same form as in the previous section (up to a real 
factor) but with /3 replaced hj 2 — p. The corresponding monic orthogonal polynomials 
Qn on the shifted lattice N + 1 — /3 satisfy 

oo 

J2Qn{k+l- P)QUk + l- P)vk = 0, n^m, (2.21) 

and they are simply the polynomials Pn{-', a,2 — f3) but shifted over a distance 1 — (3: 

Q„(x) = P„(x + /3-l;a,2-/3). (2.22) 

These are orthogonal polynomials with a positive measure whenever /3 < 2. The remark- 
able fact is that the recurrence coefficients in the three-term recurrence relation 

xQn{x) = Qn+lix) + hnQn{x) + al^Qn-l{x) (2.23) 

satisfy the same system of non-linear recurrence relations (discrete Painleve equations) as 
in Theorem 12.11 but with a different initial condition. 

Theorem 2.2. The recurrence coefficients for the orthogonal polynomials defined by 
h2.21\) satisfy the system of equations 

an 



al 



bn + bn^i -n + (3 
{al^i- a){al- a) = a(6„ - n)(6„ - n + /5 - 1), 



with initial conditions 



^I_p{2^) 
h-p{2^) ' 



where ly is the modified Bessel function 



k=0 



k\T{k + i^ + l)' 



Proof. We will give two ways to prove the result. The first way is to use the relation (12. 22 p . 
which readily gives a^ = and 6„ = 6„ + 1 — /3, where a„ and 6„ are the recurrence 
coefficients for the polynomials P„(x;a,2 — /3). We can therefore use Theorem 12. H but 
with /3 replaced by 2 — /3, and substitute by and 6„ by 6^ — 1 + /3- This indeed leaves 
the equations unchanged. The initial conditions = and 

, _ ^h-p{2^) + (1 - /3)/i-,3(2v/^) 

The latter expression can be simplified by using a well-known recurrence relation for 
modified Bessel functions [TJ Eq. 10.29.1 on p. 251]. 

An alternative way is to repeat the proof of Theorem 12. H but now on the lattice 
N + 1 — /3. Observe that we only used the Pearson equation (12. 2p and the boundary 
condition u'(— 1) = in the proof of Theorem 12.11 to arrive at the discrete equations 
|)-f l2.9p . In the present case the Pearson equation is still valid and we now use the 
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boundary condition w{—f3) = 0, which allows us to use summation by parts without the 
boundary terms. The only difference is that the moments are now 



mo 



m 



oo ^ 

oo 

1 = + 1 - + 1 - /3) 



k=0 

oo ^ 

= r(,3)a'-'E^r(^^ = r(«(v/S)-/_,(2v/S). 
SO that bo = mi/rho = ^/aI-|3{2^/a)/ Ii-j3{2y/a). □ 
2.3 Combining both lattices 

Now we can combine the two lattices and study orthogonal polynomials on the bi-lattice 
N U (N + 1 — We use the orthogonality measure /i = Cifii + €2^2, where Ci, C2 > 0, /xi 
is the discrete measure on N with weights Wk = w{k), and ^2 is the discrete measure on 
1^+1-/3 with weights Vk = w{k + l — (3). This discrete measure depends on two parameters 
Ci,C2, but the orthogonal polynomials will only depend on their ratio t = C2/C1 > 0. 
Let < /3 < 2, then /i is a positive measure and the monic orthogonal polynomials 
Rn{x) = Rn{x] a,(3,t) satisfy the three-term recurrence relation 

xRn{x) = Rn+l{x) + KRnix) + dlRn-l{x). (2.24) 

It is remarkable that these recurrence coefficients again satisfy the same non-linear re- 
currence relations (discrete Painleve equations) as in Theorem 12.11 and Theorem 12. 2^ but 
with an initial condition depending on the parameter t = C2/C1. 

Theorem 2.3. The recurrence coefficients for the orthogonal polynomials defined by 

00 00 

CiJ2f^n{k)RUk)Wk + C2j2Rn{k + l-P)RUk + l-P)Vk = 0, TTl ^ u, (2.25) 

fc=0 fc=0 

satisfy the system of equations 

bn + bn-i -n + P 



an 

~9 ' 



{al^^- a){al~ a) = a(6„ - n)(6„ - n + /3 - 1), 
with initial conditions 



al = 0, bo = a/q 



I(,i2^)+tl_^i2^) 
Ip_^{2^)+th_p{2^y 



where ly is the modified Bessel function 
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Proof. Going through the proof of Theorem 12. II we observe that only the Pearson equation 
(I2.2p is used, together with summation by parts. The boundary conditions w{—l) = 
and w{—P) = ensure that this summation by parts does not leave any boundary terms 
to evaluate at —1 (for the first sum) or — /3 (for the second sum). Hence the recurrence 
coefficients will satisfy the same non-linear recurrence relations as in Theorem 12.11 and 
12.21 The only difference is that the initial condition for bo = rhi/rho now is in terms of 
the moments 

mo = cimo + C2mo = ^{l3){Va'Y~^ {cil^-i{2y/a) + 02/1-/3(2^/0)) , 

and 

mi = cimo + C2mo = T{l3){y/a)'^~'^ [cilfs{2y/a) + C2/-/3(2A/a)) . 

□ 

We have now identified special solutions of the discrete system fl2.8p -( pl9|) with initial 
value = and 

' ^ /9_,(2Vi)+(/i-s(2yH) 

which depends on one parameter t G (0, 00). If we use the relation [TJ Eq. 10.27.2 on 
p. 251] 

2 

l-y{z) = I^{z) + - smu7iK^{z), 
n 

where is the other modified Bessel function, then the initial condition can also be 
written as 

I,3i2^) + sKf3{2^/^) 2t sin/37r 



Observe that 



Vl(2v/^) -sir;3-l(2V^)' l + t TT 



ip-i{2^) h-p{2^y 



Vi(2y^) h-p{2^y 



so that 60 (^) is a monotonically increasing function of t G [0, 00] when < /3 < 1 and a 
monotonically decreasing function when 1 < /3 < 2. One can use the Wronskian formula 
[71 Eq. 10.28.1 on p. 251] and the product formula p] Eq. 10.32.15 on p. 253] to verify 
this. For each initial value in [60 (0), 60(00)] (when < /3 < 1) or [6o(c)o), 6o(0)] (when 
1 < /9 < 2) the solution therefore corresponds to recurrence coefficients of orthogonal 
polynomials with a positive measure on the real line (in fact on the bi-lattice NU(N+1— /?)) 
whenever a > 0, and hence this solution satisfies certain positivity constraints: > for 
n >1 and 6„ > min(0, 1 — /3) for all > 0. Moreover, since the orthogonal polynomials 
are on a discrete set, we can use the familiar fact for orthogonal polynomials that between 
two zeros there is at least one point of the bi-lattice [3j Theorem 4.1 on p. 59]. The sum 
of the zeros Xi^n < X2,n < ■ ■ ■ < Xn^n of Rn hcucc Satisfies 

n n—l 

^ Xk,n >^yk 

k=l k=0 
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where |/o < Z/i < I/2 < 



are the points in the bi-lattice N U (N + 1 — i.e., 



y2k = k, y2k+i = k + 1- l3, ifO</3<l, 
y2k = k + l- /3, y2k+i = k, if 1 < /3 < 2, 
yk = k, if /3 = 1. 

The sum of the zeros is the trace of the truncated Jacobi matrix containing the recurrence 
coefficients, hence 

n n—1 
k=l k=0 

SO that we get the constraint 

2n-l 

> n{n-l)+n{l-P), if < /3 < 2 and /3 ^ 1, 

fc=0 

n-l 

> n(n-l)/2, if /3 = 1. 

k=0 

The sum Y12=o behaves hke m? /A+0{n) as n — )■ 00 when /3 7^ 1 and as 11? /2+0{n) when 
(5 = 1. This result follows from [G] Theorem 2.2] (see in particular the example of Charlier 
polynomials on p. 200). This difference in the behavior of the recurrence coefficients is 
easily explained since for /3 = 1 the two lattices coincide and we just have one lattice N. 
A similar situation appears when t — (which gives the orthogonal polynomials P„ on 
the lattice N) and when t — )■ 00, which gives the polynomials Qn on the lattice N + 1 — /3. 
We believe that these extreme cases (t = 0, t = 00 and /3 = 1) correspond to the only 
solution of with = for which > for all n > and 6„ = n + C(l) as 

n — )■ 00. 

Another special case occurs when /3 = 1/2. In this case the bi-lattice is equally spaced 
and equal to |N. The lattice points are yk = k/2 and the weight w at these points is 

w{k/2) ^ 



r((/e + l)/2)r(A;/2 + l) r(A; + 1) 

where we used Legendre's duplication formula [71 Eq. 5.5.5 on p. 138] for the gamma 
function. This means that if we take Ci = C2, then the weights of the measure /x = 
111 + /i2 are precisely the Charlier weights (Poisson distribution) with parameter 2^/a. 
The orthogonal polynomials are therefore Charlier polynomials [31 Chapter VI, §1] but 
with a scaling: 

Rn{x- a, 1/2, 1) = 2-"C'„(2x; 2^), 

where Cn{x\ a) are the monic Charlier polynomials with parameter a. The recurrence 
coefficients of Charlier polynomials Cn{x; a) are known to be = na and 6„ = n + a, so 
for the scaled polynomials 2~"C'„(2a;; a) they are = na/i and 6„ = (n + a)/2, and for 
our polynomials Rn{x) = 2~"C*„(2a;; 2y/a) we thus have 



_9 n^Ja ~ n r- ,^ 
< = ^^ bn = - + V^. (2.26) 



13 



One can indeed verify that fl2.26p is a solution of fl2.8p - fl2.9p with initial condition 

: /i/2(2v/^) + /-i/2(2v/^) _ 

/_l/2(2v/^) + /l/2(2v/^) " 



a. 



So there are special values of (3 for which the non-linear equations have simple solutions 
(as a function of n). 



3 A generalized Meixner weight 
3.1 The lattice N 

We consider the sequence {pn)n of polynomials, orthonormal with respect to the weights 

'^^ = T^^ fceN = {0,1,2,3,...}. (3.1) 
These polynomials satisfy 

oo 

where Sn,m is the Kronecker delta. These weights are positive whenever /3,7,a > 0. 
Notice that the special case (3 = 1 was studied in [1], where the authors showed that the 
recurrence coefficients a„, 6„ satisfy a limiting case of the asymmetric discrete Painleve 
equation a-dPiy. The special case /3 = 7 gives the well-known Charlier weight. We can 
write these weights in terms of a weight function w given by 

''^'^^-r(7)r(/3 + x)r(x + i) ^^'^^ 

so that Wk = w{k) for e M. 

We will prove the following result on the recurrence coefficients of the sequence of 
orthonormal polynomials. 

Theorem 3.1. The orthonormal polynomials with respect to the weights ( 13. ip on the 

lattice N, with 7,/3,a > 0, satisfy the three-term recurrence relation 

xpn{x) = a„+ip„+i(x) + hnPn{x) + a„p„_i(x), 

where the recurrence coefficients are given by the initial conditions 

_ 7aM(7+M+M0 
ao - u. Do - -5 777 — a — ^ ' 

where M{a, b, z) is the confluent hypergeometric function 

k=o ^^'f''^- 



14 



and a'^ = na — {■y — 6„ = n + 7 — /3 + a— (7 — l)t>„/a, where ('U„,fn)ngN satisfy the 
system of non-linear equations 

7-1 / 7 - /3\ 

{Un + Vn){Un+l+Vn) = — ^^^(Wn - tt) ( i;„ - a— — y 1 , (3.3) 
{Un + Vn)iUn + Vn-l) = an i^n + a) ( Un + ) ■ (3.4) 

"^^ -,,-1 V 7 J- / 

Proof. We will use the technique of ladder operators, proposed by Chen and Ismail [2]. 
The ladder operators are defined by 

A ^ \ //,^ //, ^.U{X + 1) - u{i) 
An{x) = an\pnie)Pnii~l)- ^ 7^^^, 

7^0 x+l-i 
Bn{x) = a„Vp„(£)p„_i(^-l)^^^^±ii^^u;(£), 

where u{x) = — 1 + These ladder operators satisfy 

An{x)pn-l{x) - Bn{x)pn{x) = Pn{x + 1) - Pn{x). (3.5) 

This can be seen by simplifying the left hand side using the definitions for y4„, Bn and u, 
and the Christoffel-Darboux identity 

i=o y 

After doing this we obtain 

n—l oo 

An{x)pn-lix) -Bn{x)pn{x) = Pj{x) Pr^i + l)Pj{i)w{i) . 

j=0 £=0 

Exactly the same expression is found by writing p„(x + 1) — Pn{x) as a Fourier series: 

n-l 

Pn{x + 1) - Pn{x) = an,jPj{x). 

j=0 

Equation (13. 5p gives rise to the compatibility relations 

Bn{x) + Bn+i{x) = ^— ^A„(x) - n(a; + 1) + V (3.6) 



a. 



and 



2 ^n-l(^) 



a„+iv4„+i(x) - a„ = (x - 6„)5„+i(x) - (x + 1 - bn)Bn{x) + 1. (3.7) 

An easy calculation gives that 

u{x + l)~u{e)_ 1 (7-l)(x + /3) 



x + a(7 + £-l)V 7 + x 
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Inserting this into the ladder operators, we find that 



and 



where 



a a X + 7 



Bn{x) = -r„, + -^t^, 3.9 

a ax + J 



Rn = J^PnWPnl^ - 1) . f _ M ^) 



i=0 

oo 



7 + ^-1 
7-1 



U 7+^-1 



oo 



Tn = a^VpnWPn-ll^- 1) Tw{i), 

^-^ 7 + t — 1 

£=0 ' 
oo ^ 

7 + t — 1 

We use the expressions f l3.8l) -( l3l9l) in the compatibihty relations, and compare the coeffi- 
cients of x^, X and 1. In this way we obtain a system of six equations: 

= Rn + Tn-l, (3.10) 

n n 

rn + tn+ Tn+l + U+l = iRn " KRn " KTn + /^^n + A - /? - 1 + ^ i?j + ^ , (3.11) 

3=0 j=0 
n n 

7^n + Ptn + jrn+l + = -bnlRn - bnf^Tn + 7a - /3 + 7 ^ i?,- + /3 ^ , (3. 12) 

i=o j=o 

= r„+i + - r„ - (3. 13) 

= (7 - &n)'"n+l + (/3 - &n)in+l " (l " + 7)'"n - (l - 6n + /3)in + (3.14) 



7a„+i-Rn+i + /5ci„+iT„+i — 7a„-R„_i — /3a„T„_i 

= -7&n'^n+l - PKtn+1 " 7(1 " " /^(l " &n>n + 7«- (3-15) 

By (13.101) we can substitute Rn by 1 — T^, and by (13.131) (and since ro = to = by 
definition) we can substitute r„ by —tn- If we apply these substitutions, (13. lip gives an 
expression for 6„ function of T„: 

&n = 7-(7-^)7^n + a + n-/3, (3.16) 

and (I3.12P gives 

n 

(7 - P){tn + Wi) = ihn - (7 - /3)&nT„ - 7a + /3 - 7(n + 1) + (7 - /?) J^T,. 

i=o 
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Substituting the newly found expression for bn in this equation, we obtain 

n-l 

tn + tn+1 = (7 - f3)T^ - (27 + a + n - /3 - l)r„ + 7 - 1 + ^T,-. (3.17) 

j=0 

On the other hand, ehmination of Rn and in f l3.15p gives 

7(«n+i-«n)-(7-/3)«'+iT„+i + (7-/3)a^T„_i = {^-p)bJ^^,+{^-P){l-b^)t^+^a (3.18) 
and in f l3.14p it gives 

al^^ - al = - I3){tn - tn+i) + a, (3.19) 

which after taking a telescoping sum gives us an expression for the recurrence coefficients 
Qn as a function of t„: 

al = na-{^- (3)t^. (3.20) 

Notice that since a^i > for all n, the knowledge of t„ will imply the knowledge of a„. 
Using ( I3.19P and the expression for 6„ in fl3.18p . we obtain 

a^T„__i - al^^Tn+i = (6„ + 7)Wi + (1 - &n - 7)^n, (3.21) 
in which substitution of f l3.20p gives us another relation between terms of the sequences 

{tn)n and {Tn)n- 

-a{n + l)Tn+i + anTn-1 = t„+i(n + a + 27 - /3) - t„(r;, + a + 27 - /3 - 1) 

-t„+i(7 - /3)(T„ + T„+i) + t„(7 - /3)(T„ + T„_i). 

Taking a telescopic sum, we get 

n-l 

(T„ + Tn-i){an - (7 - /3)t„) = -(a + 27 + n - /3 - l)t„ + a ^ T,-. (3.22) 

j=0 

On the other hand, if we multiply fl3.2ip by T„ and we use fl3.16p and f l3.17p . we obtain 

n n—1 

O'l+lTn+lTn - alTnTn-l = t^+i " ~ (7 " l)(^n+l - ^n) " tn+1 ^ + tn ^ 7"^, 

j=0 i=0 

which by taking a telescoping sum gives 

If we multiply fl3.22p by T„ and f l3.17p by t„, comparison of these two expressions gives 

n-l 
j=0 

Using in this equality the expression for X]j=o '^J ^^^om fl3.17p . we obtain a relation between 
Tn, tn and tn+l- 

aTn{Tn - 1)(T„(7 - /3) - 7 + 1) = (t„ + aT„)(t„+i + aT„). (3.24) 
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Finally, multiplying fl3.22p by atn, and using once again (13.231) to eliminate the sum of Tj, 
and fl3.20p to eliminate a^, we obtain a relation between T„, T„_i and tn- 

The substitution 

7 — 1 7 — 1 

allows us to write the system of two equations in a more symmetrical form: 

7-1 / 7-/3 

{Un + Vn){Un+l+Vn) = :^Vn{Vn - a) [ Vn - a - 

\ 7 — 1 

(ttn + t^n)(Mn + t^n-l) = an (^n + «) ( + f 

which is the system fl3.3p -f p7^ given in the Theorem. The initial conditions are given by 
Uq = 0, since to = 0, and Vq = ^pii'j — (3 + a — ^) since Tq and 60 are related by f l3.16p . 
and since 60 = ttT'i/^o, where is the k-th moment of the weights {wk)ke'N- The first 
two moments are easily recognized as confluent hypergeometric functions: 

00 

mo = ^ Wfe = M(7, (3, a) 

k=0 

and 

00 

mi = Y^ = ^M(7 + 1, /3 + 1, a). 

□ 

The system f l3.3p -( l37il) can be seen as a limiting case of the asymmetric discrete 
Painleve-IV equation a-dPiv ([H -Eg on p. 296], [8, Appendix A. 3]) given by 

(Y ^v\(Y - iYn-A)iYr.-B)iYr.-C)iYr.-D) 



+ Yn){Xn + Yn-l] 



+ r — Zn){Yn — F — Zn) 

(X„ + A)(X„ + B)iX„, + + D) 



{Xn + A — Z„+i/2)(X„ — A — Zn+1/2) 

where A + B + C + D = and Zn is a polynomial of degree one in n. Indeed, let 

Xn = Un ; Yn = Vn -\ , Zn = 7 I ^ ~ 77 1 H ; 

e e 7 — 1\ 2/e 



and 



A = B = a — a , C=aH — , D = — h a , 

e e 7— 1 e €7—1 



T.2 -4a^ . 2 
F = A = - 

(7-1)6' 6' 

then letting e tend to zero gives the system (I3.3p - (l3.4p . 
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3.2 The shifted lattice N + 1 - /3 

The weight function w has zeros at the negative integers —1, —2, . . ., but also at — /3, —1 — 
/3, . . .. Therefore it makes sense to consider the weight function on the shifted lattice 
N + 1 — /3. It is easy to verify that the weight on this shifted lattice is, up to a constant 
factor, equal to the weight on the original lattice N, with different parameters. Indeed, if 
we denote 

^ r(/3) r(7 + a:)a" 
^7,/3,aW r(7)r(a; + l)r(/3 + x)' 

then 

w^,fs,a{k + 1 - /3) = a r(2-/3)r(7) ^^+^-/^'2-/3,a(fc). 
Hence the corresponding orthonormal polynomials which satisfy 

oo 

<ln{k + 1 - /3)qUk + 1 - P)w{k + 1-/3) = (3.26) 

fc=0 

are, up to a constant factor, equal to the polynomials pn, shifted in both the variable x 
and the parameters 7 and (3: 

q^n^'^ix) = CnPl^'-^'^-^'^ix + (3-1). (3.27) 

The weights {w(k + 1 — f3))ke'N are positive when a>0, 7 + I — /5>0 and 2 — (3 > 0, 
i.e., whenever a > 0, (3 < 2 and 7 > /3 — 1. Remarkably, the recurrence coefficients in the 
recurrence relation 

xqn{x) = a„+ig„+i(x) + 6„g„(a;) + a„g„_i(x) 

satisfy exactly the same system of difference equations as the polynomials Pn, only the 
initial condition for 60 changes. 

Theorem 3.2. The recurrence coejficients for the orthonormal polynomials Qn, defined 
by l{3.26\) are given by = na — {•y — !)«„, 6„ = 72 + 7 — /3 + a — (7 — l)vn/a, where 
{un,Vn)n£n Satisfy the system of equations 

7 - 1 / 7-/3 

{Un + Vn){Un+l+Vn) = 5— {)„({)„- fl) {)„- O - 

a"^ \ 7 — 1 

-N/- - N "^n 1 ~ P 

{Un + Vn){Un + Vn-l) = an (^n + ( + « , 

"^n V 7 -l- 

wzt/i initial conditions 

. n A n ,, M(7-/3 + l,l-/3,a) 
ao = 0, ^o = (l-/3) ^^^_^^^^^_^^^^ , 

where M{a, b, z) is the confluent hypergeometric function. 



19 



Proof. The proof of Theorem 13.11 can easily be adapted to this case, where the ladder 
operators are now given by 



oo 



i / X , r, -, ^\ , r, ^^ + 1) — + 1 — /3) ^ ^, 

An{x) = a„2_,9n(^+l-/9)gn(^-/3)^ -P+R —^i^+^-P), 



oo 



1) 

They satisfy 

An{x)qn-iix) - 5„(x)g„(x) = g„(x + 1) - g„(x), 

leading to exactly the same compatibility relations for and Bn as in fl3.6l) - fl3.7p . Hence 
we obtain the same difference relation for the recurrence coefficients a„ and 6„ as in the 
case of the lattice N. The only difference is the initial condition for bo, which can again 
be obtained by calculating rhi/rho. 

Alternatively, if we use the relation fl3.27p . then we see that 

= a' (7 - /3 + 1, 2 - /3, a), 6„ = 6„(7 - (3 + 1,2 - (3,a) + 1 - (3, 

where a„ and 6„ are the recurrence coefficients of the polynomials in Theorem 13.11 but 
with 7 replaced by 7 — /9 + 1 and /3 replaced by 2 — /3. This means that 

7-/3 7-/3 

Un = T^nil - /3 + 1, 2 - /3, a), Vn = T^nil - (3 + 1 , 2 - (3 , tt) 

7 — 1 7 — 1 

and if we insert this in f l3.3p -f l374|) . then we retrieve the same non-linear system of equations 
for {un,Vn)nm- The initial condition is 

bo = bo{-f-P + l,2-(3,a) + l-(3, 

which becomes 

, _ a(7 - /3 + 1)M(7 + 2 - /3, 3 - /3, g) + (2 - /3)(1 - /3)M{j + 1 - /3,2 - /3,a) 
°" (2-/3)M(7-/3 + l,2-/3,a) 

This can be simplified by using some recurrence relations for confluent hypergeometric 
functions, in particular 13.3.3 and 13.3.4 in [3 p. 325], which leads to 

? ^, M(7-/3+l,l-/3,a) 
^°"^^"^V(7-/3 + l,2-/3,a)- 

□ 



3.3 Combining both lattices 

We can now consider a combination of both lattices by considering the measure /i = 
fii+tfi2 on the bi-lattice NU(N-|-1 — where /xi is the discrete measure on M with weights 
{w{k))k£n, and /i2 is the discrete measure on N + 1 — /3 with weights {w{k + 1 — f3))km- 
In order to have two positive measures /xi and /i2 we impose the conditions 

a>0, 0</3<2, 7 > max(0,/3 - 1), tG[0, oo]. 
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The orthonormal polynomials (r„)ngN for this measure depend on t > 0. Once again the 
recurrence coefficients which appear in the recurrence relation 

xr„(x) = a„+ir„+i(x) + &„r„(x) + a„r„_i(x) 

satisfy the same system of non-linear equations. 

Theorem 3.3. The recurrence coefficients for the orthonormal polynomials rn, defined 
by 

oo oo 

rn{k)r^{k)w{k) + t rn{k + 1 - (5)r^{k + 1 - (5)w{k + 1 - /3) = 6n,m 



k=0 



k=0 



are given by a„ = na — (7 — l)un and bn = ri + 7 — /3 + a— (7 — l)vn/ a, where {un, f «) 
satisfy the system of non-linear equations 



{Un + Vn){Un+l + Vn) 
{Un + Vn){Un + Vn-l) 

with initial conditions Oq = and 
where 



7-1- 



-fnU^n - a) \v„,-a 



Ur. 



u. 



[Un + a) \Un + a 



7-/9 
7-1 

7-/9 
7-1 



~ mi + tmi 

bo = —, 

mo + trrio 



7a 



(3.28) 



mo = M(7,/3,a), mi = y M(7 + 1, /3 + 1, a), 

-0 = ^r(7)r(2-^)'^ ""'^^^ - /3 + 1, 2 - /3, a), 
^1 - ;/^ a^-^M(7 - /3 + 1, 1 - A a). 

Proof. Once again, the proof of Theorem 13. II can be adapted to this case. Now the ladder 
operators are given by 



2_^rn{i)rn{i - 1) ^^—^ w{£} 



l£=0 

00 



x + 1- 



+t + 1 - p)r4i - ^ f^'^'l':^'V'^^ -i^ +1-/5) 



BJx] 



x-£ + (3 

u{x+l) -u{£) 

> r„ £ r„_i £ - 1 — — w{£) 

^-^ X + 1 — £ 



00 



00 



^=0 



They satisfy 



x-£ + /3 

AJx)rn-i(x) - BJx)rJx) = rJx + 1) - rJx) 
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leading to exactly the same compatibility relations for An and -B„ as in fl3.6l) - fl3.7l) . Hence 
we obtain the same difference relation for the recurrence coefficients a„ and 6„ as in the 
case of the lattice N or the shifted lattice N + 1 — /3. The only difference is the initial 
condition for bo, which can again be obtained by calculating rhi/rho, where rhk = mk+trhk-, 
with TUk and rhk the fc-th moment of the measures /ii and /i2 respectively. □ 

Notice that letting t tend to or oo, we obtain the initial condition for b^ of Theorem 
3.1 and 6o of Theorem 3.2 respectively, which was to be expected. If we use the formula 
O Eq. 13.2.42 on p. 325] 

^' ^) = Fr^^^T^^(«' ^' ^) + ^^^;?^^'^'M{a -b + 1,2 -b,z) 
1 (a — + Ij 1 (a) 

then the initial condition can also be written as 

r la M(7 + 1,13 + I, a)- sV{i -(3 + l)/r(-/3)[/(7 + 1,(3 + I, a) t 
bo = -7, VJ7 7^7 n . .^ i^,. ^V777 TT-. ' ^ 



(3 M(7,/3,a)-sr(7-/3 + l)/r(l-/3)f/(7,/3,a) ' 1 + t' 

where U{a, b, z) is the second solution of the confluent hypergeometric differential equa- 
tion. 

Now we have identified a class of solutions to fl3.3p - fl3.4p . with initial conditions Sq = 
and fl3.28p . depending on a parameter t. This initial condition is a homographic func- 
tion of t. Using the Wronskian formula (TJ Eq. 13.2.33] and the contiguous relations [TJ 
Eq. 13.3.15] and 

o z 

—Mia + l,b + l,z) = il- b)iMia, b, z) - Mia, 6-1, z)) 
b 

from [7, Eqs. 13.3.3 and 13.3.4], one can prove that b^ is an increasing function of t G 
]R+ when < /3 < 1, and a decreasing function of t G when 1 < /3 < 2. Hence 
for each initial value in [6o(0), 6o(oo)] or [6o(c)o), 6o(0)], the solution to this non-linear 
system corresponds to recurrence coefficients of orthogonal polynomials on the bi-lattice 
NU (N + 1 

Exactly the same results concerning the asymptotic behavior of the sum X]fc=o ^fc hold 
as in the case of the generalized Charlier polynomials: the sum behaves as n^/4 + 0{n) 
as — 7- oo when /3 7^ 1, and as ra^/2 + 0{n) when /3 = 1 or for the limiting cases t = 
and t = 00, which correspond to the latiice N and N + 1 — /3. Unfortunately, there 
seem to be no parameter choices for which the orthogonal polynomials on the bi-lattice 
turn out to be in a well-known family and for which the recurrence coefficients a^, 6„ are 
explicitly known (as was the case for (3 = 1/2, ci = C2 in the generalized Charlier case). 
Nevertheless, some parameter choices deserve special attention. The case f3 = gives 
rise to the weight function w{x) = a^ /T{x + 1). With the choice t = this gives us the 
Charlier polynomials on the lattice N. The case (3 = l/2,t = 1 gives the equally spaced 
lattice The weight in the lattice point k/2 is then given by 

'"(^•/^)-^^, 

as can be seen by using Legendre's duplication formula [71 Eq. 5.5.5]. 
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4 Concluding remarks 



Even though the recurrence coefficients of the generahzed Charher polynomials satisfy the 
same non-linear system of recurrence equations, their behavior for the bi-lattice is quite 
different from the behavior on the lattices. In Figure [1] we have plotted on the left the 
recurrence coefficients a„ and on the right the 6„ for three cases. For the a„ (left plot) 
the lowest curve corresponds to the lattice N + 1 — /3, the curve just above to the lattice 
N and the top curve corresponds to the bi-lattice with t = 10 (the parameters are a = 3, 
/3 = 1/3). For the 6„ (right plot) the top curve corresponds to the lattice N + 1 — /3, the 
curve just below to the lattice N and the lowest curve to the bi-lattice with t = 10. 




Figure 1: Recurrence coefficients for generalized Charlier polynomials (a = 3, /3 = 1/3, 
t = 10) 

Observe that the behavior for each of the two lattices is quite similar and monotonic 
as n — )■ oo, but for the bi-lattice the behavior is oscillating with a different asymptotic 
behavior as n — t- oo. We conjecture that for the lattice N and N + 1 — /3 one has the 
asymptotic behavior lim^^oofln = '^5 but for the bi-lattice one has = n^fajl + 0(1), 
where the term OiV) is bounded but oscillating. This behavior is in agreement with the 
special solution fl2.26p for /3 = 1/2. For the 6„ we conjecture that 

fo for the lattice N, 
lim On — n = < 

[1-13 for the lattice N + 1 - /3, 

and bn = n/2 + 0{1) for the bi-lattice. 

Similar observations hold for the recurrence coefficients of the generalized Meixner 
polynomials. In Figure |2] we have plotted the recurrence coefficients a„ (on the left) and 
bn (on the right) for the same three cases for the parameter values a = 3, /3 = 2/3, 
7 = 9/10. For the a„, the lowest curve corresponds to the lattice N-|- 1 — /3, the curve just 
above to the lattice N (both curves are almost identical), and the top curve corresponds 
to the bi-lattice with t = 2. For the bn the top curve corresponds to the lattice N-|- 1 — /3, 
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Figure 2: Recurrence coefficients for generalized Meixner polynomials {a = 3, /3 = 2/3, 



the one just below to the lattice N (both curves are almost identical), and the lowest 
curve corresponds to the bi-lattice. Here we conjecture that for each lattice one has 



C(l) and bn/n = C(l), where the C(l) terms are oscillatory. 
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